Rf-induced persistent long-range ordered structures in two-species ion Coulomb 

crystals in a linear Paul trap. 
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We report on the observations of three-dimensional long-range ordered structures in the central 
40 Ca + ion component of 40 Ca + - 44 Ca + two-species ion Coulomb crystals in a linear Paul trap. In con- 
trast to long-range ordering previously observed in single species crystals, the structures observed 
are strikingly persistent and always of one specific type in one particular orientation. Molecular 
dynamics simulations strongly indicate that these characteristics are a hitherto unpredicted conse- 
quence of the co-axial cylindrical symmetry of the central ion component of the Coulomb crystal 
and the radio frequency quadrupole trapping field. 

PACS numbers: 32.80.Pj, 52.27.Jt, 52.27.Gr, 36.40.Ei 
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A solid state of an one-component plasma (OCP), 
often referred to as a Wigner crystal or a Coulomb 
crystal, appears whenever the coupling parameter T = 
^Coui /-Ekin > where i?coui is the nearest neighboring 
Coulomb potential energy and £"kin is the averaged ki- 
netic energy of the particles, exceeds ~200 QHH. In re- 
cent decades, such crystals have been investigated exper- 
imentally in a large variety of physical systems [^IfjIaLj]. 
In particular, laser-cooled and trapped atomic ions have 
proven to be an excellent system for experimental stud- 
ies of Coulomb crystals under various confinement con- 
ditions. Single-species, one-, two- and three-dimensional 
ion Coulomb crystals have been investigated in Penning 
and radio frequency (rf) traps (also named Paul traps), 
with structural findings in good agreement with theo- 
retical predictions 0, H H [3 13 111 El ■ Two-species 
crystals have been studied in much less detail, and exper- 
imentally mainly in linear Paul traps [lii . Il4| . For singly 
charged ions in such traps, the lighter species will in gen- 
eral segregate into a cylindrically shaped structure sur- 
rounded by the heavier ion species. In the case of 40 Ca + 
and 24 Mg + ions, the structures of the lighter species ^3] 
were found for the most part to be identical to the cylin- 
drical structures of infinitely long one-species systems 
confined in two-dimensions by a rotational-symmetric 
harmonic potential [To] ]. For two ion species with iden- 
tical c harg e-to- mass ratios, mixing of the species is pre- 
dicted At present, Coulomb crystals find applica- 
tions in such diverse fields as quantum computing [l7lll8j | 
and cold molecule ion research [ldl liH l2flj . 

In this Letter, we present observations of persistent 
long-range ordered structures in the 40 Ca 4 " component of 
bicrystals consisting of 40 Ca + and 44 Ca + ions. The ob- 
served structures deviate from those expected to exist in 
a fully rotational symmetric harmonic potential [lii . Ilfij , 
but are in close agreement with results from molecular 
dynamics (MD) simulations including the quadrupole na- 
ture of the trapping fields. 

The Ca + ions are confined in a linear Paul trap which is 



situated in a vacuum chamber at a pressure of 10 -10 Torr, 
and are laser cooled as depicted in Fig. 1(a). The Paul 
trap used in these experiments has been described in de- 
tail elsewhere [2lJ. In short, the linear Paul trap consists 
of four electrode rods placed in a quadrupole configura- 
tion. The electrode diameter is 8.0 mm and the minimum 
distance of the central trap axis is r = 3.5 mm. Time 
varying voltages gf/rf cos(f2 r ft) and ^U r f cos(£l r fi + 7r) are 
applied to the two sets of diagonally opposite electrode 
rods, respectively, to obtain confinement in the radial 
plane (xy-plane in Fig. 1(a)). Axial confinement along 
the z-axis is accomplished by sectioning each of the elec- 
trode rods into three pieces and then applying a static 
voltage [/end to the end-electrodes. The length of the 
center-electrode is 5.4 mm, while the outer pieces are 20 
mm. In the present experiments fi r f = 2ir x 3.88 MHz, 
Uri ~ 540 V and U cnd ~ 10-50 V were used. The result- 
ing confinement for an ion species of type i with mass 
Mi and charge Qi is often approximated by a harmonic 
pseudo potential $ ps (r, z) = ^Mi(ui 2 r 2 + oj 2 z 2 ), where 
lo t and w z are the radial and axial trap frequencies, re- 
spectively. The axial trap frequency is given by u 2 = 
2KQiU en d/Mi, where k = 3.97 x 10 4 mT 2 is a constant 
related to trap geometry, and the radial trap frequency 
is given by uj 2 r = uj 2 f - \u 2 z , where to 2 t = Q 2 U^/2M i r%Q, 2 i 
is the contribution from the time varying quadrupole 
fields. The dependence on the charge and mass of the ion 
species makes the lighter isotope ( 40 Ca + ) more tightly 
bound towards the trap axis than the heavier ( 44 Ca + ) 
and leads consequently to a total radial separation of the 
two ion species when sufficiently cooled [ 1 3j . The zero- 
temperature ion density in the pseudopotential is given 
by nthoo = £oU 2 f /MirQ^l 2 . { , where eo is the vacuum per- 
mittivity 0] • Due to the spatial separation of the ions in 
two-species ion Coulomb crystals, this expression is also 
applicable to the individual components of such crystals. 

The 40 Ca + and 44 Ca + ions used in the experiments 
are produced isotope selectively by resonant two-photon 
photo-ionization of atoms in an effusive beam of naturally 
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FIG. 1: (a) Schematics of the laser and linear Paul trap setup. 
(D)M denotes a (dichroic) mirror, PBS denotes polarization 
beamsplitters, and DL denotes diode laser, (b) Optical tran- 
sitions used for Doppler laser cooling of the Ca + ions. 



abundant calcium |22j,l23j. In Fig. 1(b), the transitions 
in Ca + used for Doppler laser cooling of the trapped ions 
are shown. The main cooling transition is the dipole al- 
lowed 4Si/ 2 — > 4Pi/2 transition at 397 nm. To avoid opti- 
cal pumping into the metastable 3D 3 /- 2 state, rcpumping 
is done either by using a single repump laser at 866 nm 
via the 4P 1 / 2 state ( 44 Ca + ) or by using two repump lasers 
at 850 nm and 854 nm via the 4P 3/2 state ( 40 Ca+L Due 
to the isotope shifts of the cooling transitions [3, , 
each isotope ion requires its own laser cooling frequencies. 
For both isotopes the final temperature is of the order 
of ~10 mK, which is low enough for achieving Coulomb 
crystallization (r ~ 250 0]). Imaging of the fluorescence 
from the trapped ions is achieved using an image inten- 
sified CCD camera placed above the trap. 



FIG. 2: (color online) Images of two-species ion Coulomb 
crystals containing ~1500 40 Ca+ ions (red) and ~20 44 Ca+ 
ions (blue) at different settings of the end-cap potential, U cn d- 
The camera exposure time is ~100 ms and the trap potentials 
are U lt = 540 V and (a) U end = 46.1 V, (b) U md = 30.2 V, 
(c) (/end = 20.0 V, and (d) U cn d = 13.8 V, respectively. 
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Images of two-species Coulomb crystals consisting of 
40 Ca + and 44 Ca + ions are presented in Fig. 2 for various 
settings of the static voltage C/ cnd on the end-electrodes. 
Since the ions are only fluorescing when they are directly 
laser cooled, the two isotopes can be imaged separately 
by alternatingly blocking (about 1 Hz rep. rate) the 850 
nm ( 40 Ca + laser cooling off) and the 866 nm laser ( 44 Ca + 
laser cooling off) using mechanical shutters. The pre- 
sented combined images are subsequently created with 
a red color coding for 40 Ca + and blue for 44 Ca + . Due 
to sympathetic cooling [2(| [2?|], the crystal retains its 
shape and structure during the alternating laser cooling 
sequence. As expected [l3j, the lightest isotope 40 Ca + 
is located as a cylindrical core closest to the trap axis, 
surrounded by the heavier 44 Ca + . It is immediately clear 
from the images in Fig. 2 that the projection of the ac- 
tual three-dimensional structure of the 40 Ca + ions is a 
two-dimensional rectangular lattice aligned with the trap 
axis. Since the depth of focus of the imaging system (~50 
fim) is several times the inter-ion distance, we conclude 
the rectangular structure in the images must originate 
from a three-dimensional long-range ordering || . In con- 
trast to our previous observations of long-range order in 
spherical one-component ion crystals || , where the orien- 
tation of the observed metastablc (~100 ms) structures 
seemed to be arbitrary, the rectangular structures pre- 
sented here are very persistent (~10 s) and always ori- 
ented the same way. Accordingly, the presence of the sur- 
rounding 44 Ca + ions apparently has significant influence 
on the formation and appearance of the observed long 
range structure in the 40 Ca + part of the crystal. The 
image sequence in Fig. 2 illustrates additionally that de- 
spite changes in the outer shape of the 40 Ca + core, the 
observed rectangular lattice of the ions is preserved, in- 
dicating that the observed long-range ordered structure 
is rather stable to changes in the boundary conditions of 
the crystal. 

In order to understand the observations, a series of 
molecular dynamics (MD) simulations of two-species 
crystals with the same number of the two calcium iso- 
tope ions as in the crystals shown in Fig. 2 have been 
performed. In Figs. 3(a) and 3(b), results from one sim- 
ulation using a pseudo-potential corresponding to the 
trapping parameters of Fig. 2(b) is presented. As in 
the experiments, a clear radial separation of the two iso- 
tope ions is observed. Furthermore the radial projec- 
tion (Fig. 3(a)) clearly reveals that the 40 Ca + part of 
the crystal is organized in concentric cylindrical struc- 
tures, resembling the structure of an infinitely l ong ID 
cylindrically symmetric confined ion crystal |l3L lla j . In 
Fig. 3(b), a projection corresponding to the focal re- 
gion of the imaging system is shown. Neither this nor 
other sections, e.g., in the j/z-plane, lead to projection 
images with rectangular structures. However, when the 
full rf potential is used in the MD simulations, some 
regular ordering in the central component of the crys- 



FIG. 3: (color online) MD simulation data of a two-species 
ion Coulomb crystal containing 1500 40 Ca + ions (red) and 
2000 44 Ca + ions (blue) using a pseudo potential approxima- 
tion ((a) and (b)) and using the full rf potential ((c) and (d)), 
respectively. While (a) and (c) show a projection of all the 
ion positions onto the zy-plane, (b) and (d) present the data 
points contained in a slice of thickness 24 (im through the 
crystal center projected to the 22-plane. The slices are indi- 
cated by white dashed lines in (a) and (c). Trap potentials 
are U rt = 540 V and U cnd = 33 V. 



tal does appear, as is evident from the results presented 
in Figs. 3(c) and 3(d) for a specific phase of the rf field. 
From Fig. 3(d), it is seen that indeed a rectangular pro- 
jection image is expected when the rf-quadrupole field, 
which breaks the rotational symmetry, is included in the 
simulations. Even when averaging over all phases of the 
rf field, the rectangular structure persists. However, some 
blurring of the position of the ions along the x-axis, as is 
seen in the images of Fig. 2, is found. Analysis of a much 
simpler two-ion system in a linear rf trap has previously 
shown similar preferred orientation effects with respect 
to the rf quadrupole field axes |2gj . 

A closer analysis of the simulation results presented 
in Figs. 3(c) and 3(d) shows that the 40 Ca + ions in the 
core organize themselves in a long-range ordered struc- 
ture in the form of a face-centered tetragonal (fct) lat- 
tice as illustrated in Fig. 4(a), with the sidclcngths re- 
lated by a = b = \/3c. The rectangular projection pre- 
sented in Fig. 3(d) is obtained when the fct structure 
is viewed along the b vector as illustrated in Fig. 4(b). 
This rectangular projection has a height to width ratio 
of hfet/wfet = V3 ~ 1.73, which is not exactly the same 
as the h/w = 1.62 ±0.07 observed in the experiments. In 
fact, the observed rectangular structure is more in agree- 
ment with a face-centered cubic (fee) structure viewed 
along the [211] direction (see Figs. 4(c) and 4(d)), which 
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FIG. 4: Fct lattice cell (a) with projection along the a vec- 
tor (b). Fee lattice cells (c) with projection along the ^In- 
direction (d). 



would lead to h{ cc /w{ cc = •v/8/3 ~ 1.63. The reason for 
the deviation between the simulated results and the ac- 
tually observed structures is probably that the difference 
in the potential energies of the two structures in the rf 
potential is very small, as is well-known to be the case for 
various long-range ordered structures in infinite systems 
without the presence of rf fields Q . In recent single com- 
ponent experiments both body-centered cubic (bec) and 
fcc-like structures were indeed observed @, but in con- 
trast to the two-species results above, no fixed orientation 
with respect to the trap axis was found. Another point 
supporting that the observed structure is a fee structure 
is the ion density. Assuming that the observed projection 
images of 40 Ca + ions in Fig. 2 are actually fee structures 
observed along the [211] direction, the ion density must 
be nf cc = (3.8 ± 0.4) x 10 8 cm -3 . In comparison, the 
40 Ca + ion density calculated from the trap parameters is 
n thoo,40 = (4-3 ± 0.3) x 10 s , in good agreement with the 
fee assumption. 

Very stable aligned crystal structures as those dis- 
cussed above, may find many future applications. For 
instance, for cavity QED studies, the situation where 
the axial periodicity of the Coulomb crystal is an integer 
multiple of the node spacing of the standing wave cav- 
ity field is very interesting, since the effective coupling of 
the atomic ensemble to the light field can be controlled 
by shifting the position of the whole crystal. 

In conclusion, very persistent three-dimensional long- 
range ordered structures with one specific orientation 
have been observed in two-species ion Coulomb crystals 
in a linear Paul trap. MD simulations strongly indi- 
cate that these characteristics are a consequence of the 



co-axial cylindrical symmetry of the central ion compo- 
nent of the Coulomb crystal and the radio frequency 
quadrupolc trapping field. 
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